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Abstract 

This work studies the generalized Moran process, as introduced by Lieberman, Hauert, and 
Nowak [Nature, 433:312-316, 2005], where the individuals of a population reside on the vertices of an 
undirected connected graph. The initial population has a single mutant of a fitness value r (typically 
r > 1), residing at some vertex v of the graph, while every other vertex is initially occupied by an 
individual of fitness 1. At every step of this process, an individual (i.e. vertex) is randomly chosen 
for reproduction with probability proportional to its fitness, and then it places a copy of itself on 
a random neighbor, thus replacing the individual that was residing there. The main quantity of 
interest is the fixation probability, i.e. the probability that eventually the whole graph is occupied by 
descendants of the mutant. In this work we concentrate on the fixation probability when the mutant 
is initially on a specific vertex v, thus refining the older notion of Lieberman et al. which studied the 
fixation probability when the initial mutant is placed at a random vertex. We then aim at finding 
graphs that have many "strong starts" (or many "weak starts") for the mutant. Thus we introduce a 
parameterized notion of selective amplifiers (resp. selective suppressors) of evolution, i.e. graphs with at 
least some h(n) vertices (starting points of the mutant), which fixate the graph with large (resp. with 
small) probability. We prove the existence of strong selective amplifiers (i.e. for h(n) = O(n) vertices 
v the fixation probability of v is at least 1 — ^£li for a function c(r) that depends only on r). We also 
prove the existence of quite strong selective suppressors. Regarding the traditional notion of fixation 
probability from a random start, we provide the first non-trivial upper and lower bounds: first, we 
demonstrate the non-existence of "strong universal" amplifiers, i.e. we prove that for any undirected 
graph the fixation probability from a random start is at most 1 — ^374- Finally we prove the Thermal 
Theorem, which states that for any undirected graph, when the mutant starts at vertex v, the fixation 
probability at least (r — l)/(r + d ^° gl * ). This implies the first nontrivial lower bound for the usual 
notion of fixation probability, which is almost tight. This theorem extends the "Isothermal Theorem" of 
Lieberman et al. for regular graphs. Our proof techniques are original and are based on new domination 
arguments which may be of general interest in Markov Processes that are of the general birth-death type. 
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1 Introduction 

Population and evolutionary dynamics have been extensively studied [2j[6j[7j[T5|[2l][24j[25], mainly on the 
assumption that the evolving population is homogeneous, i.e. it has no spatial structure. One of the main 
models in this area is the Moran Process |l9j, where the initial population contains a single mutant with 
fitness r > 0, with all other individuals having fitness 1. At every step of this process, an individual is 
chosen for reproduction with probability proportional to its fitness. This individual then replaces a second 
individual, which is chosen uniformly at random, with a copy of itself. Such dynamics as the above have 



been extensively studied also in the context of strategic interaction in evolutionary game theory 11 14 23 



In a recent article, Lieberman, Hauert, and Nowak [16] (see also 20 ) introduced a generalization of the 
Moran process, where the individuals of the population are placed on the vertices of a connected graph 
(which is, in general, directed) such that the edges of the graph determine competitive interaction. In the 
generalized Moran process, the initial population again consists of a single mutant of fitness r, placed on a 
vertex that is chosen uniformly at random, with each other vertex occupied by a non- mutant of fitness 1. An 
individual is chosen for reproduction exactly as in the standard Moran process, but now the second individual 
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to be replaced is chosen among its neighbors in the graph uniformly at random (or according to some weights 
of the edges) (l6|[20]. If the underlying graph is the complete graph, then this process becomes the standard 



Moran process on a homogeneous population 16|20 . Several similar models describing infections and particle 



interactions have been also studied in the past, including the SIR and SIS epidemics 10 Chapter 21], the 
voter and antivoter models and the exclusion process [TJ|9j[l7]. However such models do not consider the 
issue of different fitness of the individuals. 

The central question that emerges in the generalized Moran process is how the population structure affects 



evolutionary dynamics 16 20 . In the present work we consider the generalized Moran process on arbitrary 
finite, undirected, and connected graphs. On such graphs, the generalized Moran process terminates almost 
surely, reaching either fixation of the graph (all vertices are occupied by copies of the mutant) or extinction of 
the mutants (no copy of the mutant remains). The fixation probability of a graph G for a mutant of fitness r, 
is the probability that eventually fixation is reached when the mutant is initially placed at a random vertex 
of G, and is denoted by f r {G). The fixation probability can, in principle, be determined using standard 
Markov Chain techniques. But doing so for a general graph on n vertices requires solving a liner system 
of 2™ linear equations. Such a task is not computationally feasible, even numerically. As a result of this, most 
previous work on computing fixation probabilities in the generalized Moran process was either restricted to 
graphs of small size [6] or to graph classes which have a high degree of symmetry, reducing thus the size of 
the corresponding linear system (e.g. paths, cycles, stars, and cliques [3| |5]) . Experimental results on the 
fixation probability of random graphs derived from grids can be found in [22] . 

A very interesting recent result [8] shows how to construct fully polynomial randomized approximation 
schemes (FPRAS) for the probability of reaching fixation (when r > 1) or extinction (for all r > 0). The 
result of [8] uses a Monte Carlo estimator, i.e. it runs the generalized Moran process several timea^J while each 
run terminates in polynomial time with high probability [8]. Note that improved lower and upper bounds 
on the fixation probability immediately lead to a better estimator here. Ontil now, the only known general 
bounds for the fixation probability on connected undirected graphs, are that f r (G) > ~ and f r (G) < 1— 
these are quite weak bounds and easy to derive. 

Lieberman et al. [16[|20] proved the Isothermal Theorem, namely that the fixation probability of a "sym- 
metric directed" graph is equal to that of the complete graph (i.e. the homogeneous population of the 
standard Moran process), which tends to 1 — £ when the size n of the population grows. Intuitively, in the 
Isothermal Theorem, every vertex of the graph has a temperature which determines how often this vertex 



is being replaced by other individuals during the generalized Moran process. In the terminology of 16 
symmetric directed graphs become regular graphs (i.e. graphs with overall the same vertex degree) when we 
consider undirected graphs. Furthermore, the complete graph (or equivalently, any regular graph) serves as 
a benchmark for measuring the fixation probability of an arbitrary graph G: if f r (G) is larger (resp. smaller) 
than that of the complete graph then G is called an amplifier (resp. a suppressor) |16||20| . Until now 
only graphs with similar (i.e. a little larger or smaller) fixation probability than regular graphs have been 
identified [3j-[5| [T6|[T8] , while no class of strong amplifiers/suppressors is known so far. 

Our contribution. The structure of the graph, on which the population resides, plays a crucial role in the 
course of evolutionary dynamics. Human societies or social networks are never homogeneous, while certain 
individuals in central positions may be more influential than others (20j. Motivated by this, we introduce 
in this paper a new notion of measuring the success of an advantageous mutant in a structured population, 
by counting the number of initial placements of the mutant in a graph that guarantee fixation of the graph 
with large probability. This provides a refinement of the notion of fixation probability in undirected graphs. 
Specifically, we do not any more consider the fixation probability as the probability of reaching fixation when 
the mutant is placed at a random vertex, but we rather consider the probability f r (v) of reaching fixation 
when a mutant with fitness r > 1 is introduced at a specific vertex v of the graph; f r (v) is termed the 
fixation probability of vertex v. Using this notion, the fixation probability f r (G) of a graph G = (V, E) with 
n vertices is f r (G) = £ J2vev fr( v )- 

We aim in finding graphs that have many "strong starts" (or many "weak starts") of the mutant. Thus 



1 For approximating the probability to reach fixation (resp. extinction), one needs a number of runs which is about the inverse 
of the best known lower (resp. upper) bound of the fixation probability. 
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we introduce the notions of (h(n), g(n)) -selective amplifiers (resp. (h(n), g(n))-selective suppressors), which 
include those graphs with n vertices for which there exist at least h(n) vertices v with f r (v) > 1 — 

(resp. f r (v) < ^fer) f° r an appropriate function c(r) of r. We contrast this new notion of (h(n), g(n))- 
selective amplifiers (resp. suppressors) with the notion of g(n)- universal amplifiers (resp. suppressors) which 
include those graphs G with n vertices for which f r (G) > 1 — (resp. f r (G) < jjds) for an appropriate 
function c(r) of r. For a detailed presentation and a rigorous definition of these notions we refer to Section [2] 

Using these new notions, we prove that there exist strong selective amplifiers, namely (O(n), rt)-selective 
amplifiers (called the urchin graphs). Furthermore we prove that there exist also quite strong selective 
suppressors, namely ( > ^-^jy )-selective suppressors (called the (f>(n) -urchin graphs) for any function 
4>{n) — with 4>(n) < i/n. 

Regarding the traditional measure of the fixation probability f r (G) of undirected graphs G, we provide 
the first non-trivial upper and lower bounds, while only the trivial bounds i and 1 — were known J8J. 
More specifically, first of all we demonstrate the nonexistence of "strong" universal amplifiers by showing 

that for any graph G with n vertices, the fixation probability f r (G) is strictly less than 1 3/4+s 1 f° r an y 

£ > 0. This is in a wide contrast with what happens in directed graphs, as Lieberman et al. |16| provided 
directed graphs with arbitrarily large fixation probability (see also |20|). 

On the other hand, we provide our lower bound in the Thermal Theorem, which states that for any vertex 
v of an arbitrary undirected graph G, the fixation probability f r (v) of v is at least (r—l)/(r+ d ^° s " ) for any 
r > 1, where degv is the degree of v in G (i.e. the number of its neighbors) and deg min (resp. deg max ) is the 
minimum (resp. maximum) degre in G. This result extends the Isothermal Theorem for regular graphs [16] . 
In particular, we consider here a different notion of temperature for a vertex than [16] : the temperature of 
vertex v is g^j- As it turns out, a "hot" vertex (i.e. with hight temperature) affects more often its neighbors 
than a "cold" vertex (with low temperature). The Thermal Theorem, which takes into account the vertex v 
on which the mutant is introduced, provides immediately the first non-trivial lower bound (r — 1) /(r+ ^ max ) 
for the fixation probability f r (G) of any undirected graph G. The latter lower bound is almost tight, as 
it implies that f r (G) > for a regular graph G, while the Isothermal Theorem implies that the fixation 
probability of a regular graph G tends to as the size of G increases. Note that our new upper/lower 
bounds for the fixation probability lead to better time complexity of the FPRAS proposed in j8j, as the 
Monte Carlo technique proposed in [8] now needs to simulate the Moran process a less number of times (to 
estimate fixation or extinction). 

Our techniques are original and of a constructive combinatorics flavor. For the class of strong selective 
amplifiers (the urchin graphs) we introduce a novel decomposition of the Markov chain M of the generalized 
Moran process into n — 1 smaller chains Aii, M.2, ■ ■ ■ , M. n -i, and then we decompose each Aik into two 
even smaller chains M\,M."j,. Then we exploit a new way of composing these smaller chains (and returning 
to the original one) that is carefully done to maintain the needed domination properties. For the proof 
of the lower bound in the Thermal Theorem, we first introduce a new and simpler weighted process that 
bounds fixation probability from below (the generalized Moran process is a special case of this new process) . 
Then we add appropriate dummy states to its (exponentially large) Markov chain, and finally we iteratively 
modify the resulting chain by maintaining the needed monotonicity properties. Eventually this results to the 
desired lower bound of the Thermal Theorem. Finally, our proof for the non-existence of strong universal 
amplifiers is done by contradiction, partitioning appropriately the vertex set of the graph and discovering 
an appropriate independent set that leads to the contradiction. 

The paper is organized as follows. Preliminaries and notation are given in Section [2] Furthermore we 
present our results on amplifiers and suppressors in Sections [3] and [4] respectively. 

2 Preliminaries 

Throughout the paper we consider only finite, connected, undirected graphs G = (V,E). Our results apply 
to connected graphs as, otherwise, the fixation probability is necessarily zero. The edge e £ E between two 
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vertices u, v <E V is denoted by e = uv. For a vertex subset X C V , we write X + z/ and I-yforlU {y} 
and Xn{j/}, respectively. Furthermore, throughout r denotes the fitness of the mutant, while the value r is 
considered to be independent of the size n of the network, i.e. we assume that r is constant. For simplicity 
of presentation, we call a vertex v "infected" if a copy of the mutant is placed on v. For every vertex subset 
S CV we denote by f r (S) the fixation probability of the set S, i.e. the probability that, starting with exactly 
jS*! copies of the mutant placed on the vertices of S, the generalized Moran process will eventually reach 
fixation. By the definition of the generalized Moran process / r (0) = and f r (V) = 1, while for S ^ {0, V}, 

, _ Z^efi {sk^JiS + y) + jj^fjS - a)) 

V I — £ I i— ) 

Z-~ixy£E I dcg a; dcg y J 

Therefore, eliminating self-loops in the above Markov process, 

. , _ Sxyeg^xGS^gS (dc^M-S' + T/) + dej^/^ 51 -^)) 
Exye-E,xe5,i/^5 (^degx ^~ dcg y 

In the next definition we introduce the notions of universal and selective amplifiers. 

Definition 1 Let Q be an infinite class of undirected graphs. If there exists an no G N, an rg > 1, and some 
function c(r), such that for every graph G £ Q with n > tiq vertices and for every r > ro-' 

• fr{G) > 1 — j^j? then Q is a class of g(n)-universal amplifiers, 

• there exists a subset S of at least h(n) vertices of G, such that f r (v) > 1 — for every vertex v G S, 
then Q is a class of {h(n), <7(n))-selective amplifiers. 



Moreover, Q is a class of strong universal ( resp. strong selective ) amplifiers if Q is a class of n-universal 
(resp. (0(r7,), n)-selective / ) amplifiers. 

Similarly to Definition [T] we introduce the notions of universal and selective suppressors. 

Definition 2 Let Q be an infinite class of undirected graphs. If there exist functions c{r) and no(r), such 
that for every r > I and for every graph G € Q with n > no(r) vertices: 

• fr(G) < then Q is a class of (jr(n)-universal suppressors, 

• there exists a subset S of at least h(n) vertices of G, such that f r (v) < for every vertex v G S, 
then Q is a class of (h(n), g(n))-selective suppressors. 



Moreover, Q is a class of strong universal (resp. strong selective,) suppressors ifQ is a class of n-universal 
(resp. (0(n), n)-selectivej suppressors. 

Note that ?i = n (r) in Definition [2j while in Definition [I] n is not a function of r. The reason for this 
is that, since we consider the fitness value r to be constant, the size n of G needs to be sufficiently large with 
respect to r in order for G to act as a suppressor. Indeed, if we let r grow arbitrarily, e.g. if r = n 2 , then for 
any graph G with n vertices the fixation probability f r {v) tends to 1 as n grows. The next lemma follows 
by Definitions [l] and [2j 

Lemma 1 IfQ is a class of g(n)-universal amplifiers (resp. suppressors), then Q is a class of (O(n), g{n))- 
selective amplifiers (resp. suppressors). 
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Proof. Suppose that Q is a class of g(n)-universal amplifiers. That is, for every r > tq and for every 

c ( r ) 

g( n ) 



graph G = (V, E) £ Q with n > n vertices, the fixation probability of G is f r (G) > 1 — ^f\, where c(r) 



is some function that depends only on r. Let S C V be the subset of vertices such that f r (v) > 1 — 

for some function c'(r) that depends only on r. Then there exists an appropriate function <j)(n,r) = w(l), 

i.e. lim <j)(n, r) — oo, such that f r (v) < 1 — for every u£ 7\S. Thus the fixation probability of G is 

^ J^I^ + ^-l^l)^ 1 -^) ,, (»-|g|) fln.r) 
n n g( n ) 

Now, since / r (G) > 1 - fg, it follows by (B that (n - |5|) < n^gy, and thus \S\ > n(l - ^gy) = 0(n), 
since (f>(n,r) = uj(1). Thus it follows by definition of the set 5 that Q is a class of (0(n), g(n))-selective 
amplifiers. 

Suppose now that Q is a class of <7(n)-universal suppressors. That is, for every r > 1 and for every graph 
G = (V,E) € Q with rt > no(r) vertices, the fixation probability of G is f r (G) < ^fejj where c(r) is some 

function that depends only on r. Let S C V be the subset of vertices such that / r (f) < for some function 
c'(r) that depends only on r. Then there exists an appropriate function </>(n, r) = w(l), i.e. lim <j)(n, r) = oo, 

0(n,r) 
9(n) 



such that /r(w) > ^j^y f° r every w £ 7 \ 5. Thus the fixation probability of G is 



/ (g) > I^O+^H^Jg 1 = (n-jgj) 0(n,r) 
n n g( n ) 

Now, since / r (G) < fg, it follows by (g) that (n - |5|) < n^gy, and thus |5| > n(l - ^gy) = 0(n), 
since <p(n,r) — w(l). Thus it follows by definition of the set 5 that £/ is a class of (0(n), g(rt))-selective 
suppressors. ■ 

The most natural question that arises by Definitions [l] and [2] is whether there exists any class of strong 
selective amplifiers/suppressors, as well as for which functions h{n) and g(n) there exist classes of g(n)- 
universal amplifiers/suppressors and classes of (h(n), g(n))-selective amplifiers/suppressors. In Section [3] 
and [4] we provide our results on amplifiers and suppressors, respectively. 



3 Amplifier bounds 



In this section we prove that there exist no strong universal amplifiers (Section 3.1), although there exists a 
class of strong selective amplifiers (Section pO]). 



3.1 Non-existence of strong universal amplifiers 

Theorem 1 For any function g(ri) = £l(ni +£ ) for some e > 0, there exists no graph class Q of g(n) -universal 
amplifiers for any r > r*o = 1 . 

Proof. Let tq = 1 and g(n) — i7(n 1_(5 ), where S = \ — e < 4. Suppose that Q is a class of <?(n)-universal 
amplifiers. That is, for every graph G = (V, E) € Q with n > tiq vertices, the fixation probability of G is 
fr{G) > 1 — ^g > 1 — for every r > 1, where c(r), cq(t') are two functions that depend only on r. We 
partition the vertex set V into two three sets Vx , Vi , V3 such that 

Vx = {veV:f r (v)>l-^} (4) 

V 2 = { veV \Vx:fr(v)>l-^} (5) 

V 3 = {veV\(VxUV 2 ):f r (v)<l-^±} (6) 
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where ci(r) is an appropriate function of r, and lim <p(n,r) — oo, i.e. <p(n,r) = w(l). Note that Vi ^ 0, 

n— >oo 

since f r (G) > 1 — ^rj- by assumption. Using (jij), the fixation probability f r (G) of G is upper bounded by 



/r(G) < 



\V 3 \ <j>(n,r) 



Now, since f r (G) > 1 - it follows by |7| that 1 - ^ < l 



,1-25 



|Vbl 0fer) 



and thus 



IVs I < n 



1-5 °o(r) 



(J) 



(8) 



For an arbitrary vertex v G V, we obtain an upper bound of the probability f r (v) by assuming that 
fixation is reached if the process reaches at least two infected vertices, when it starts with only v being 
infected. Therefore 



fr{v) < 



r ' 1 + J2x£N(v) degx '® 



r + EieATfij) dog a; r + Sa:£JV('u) deg x 

for every d6F. It follows now by Q and Q that for every u € V%, 



(9) 



1 - 4^4 < 1 



£N(v) deg a; 



,1-5 



y - 



xeN(v) 



degx 



< 



r + 12x£N(v) degx 

r • co(r) < c'(r) 



c (r) ft 



i-<5 



for an appropriate function c'(r) of r. Therefore, since 2xeAT(V) degx — deg(i>) ' «' (l^M implies that 

deg v < c'(r) ■ n s 

for every v G V x . Furthermore, since < J2 x eN(v) dej^ for ever y u £ N ( v )> (10) im P lies that 

degw > — T 
c(r) 

for every u G N(v), where u G Vi. Similarly it follows by ^ and (|9| that for every v G V2, 

y — 



for some function c"(r) of r, and that 



< r-ci(r) < c"(r) 
deg it ~~ n 1_2<5 — ci(r) ~~ n 1 ^ 28 



deg i; < c"(r) ■ n 

n l-2<5 

deg u > —— 
c'(r) 



26 



(10) 

(11) 

(12) 
(13) 

(14) 
(15) 



for every v G V2 and u G N{v). 

Consider now a vertex v G V± and a vertex u G N(v). Then, since 6 < \ by assumption, (11), (12 1, 
and ( 14) imply that u ^ V\ and u £ V%, and thus u G V3. That is, iV(u) C V3 for every v G Vi, an<TTi< is an 
independent set. 



G 



Let now v € V such that f r (v) is maximized. Then clearly v € V±. Denote Q v = 
it 

^2x€N(u)\{v} dcg(x) ■ 



jxEN(v) dcg(x) ' 

that Q v = n(^). Furthermore, denote for any neighbor u G N(v) the quantity Q uv = J2xeN(v)\{u} deg(x) 
1 Then it follows by 



note 



18 



that 



f r (G) < max 



2r 2 



ueN(v) 2r 2 + Q v Quv 



(16) 



cf. Theorem 1 in 18 and its proof therein. Let uq € N(v) be such that the right hand side of (16 1 is 
maximized, and thus 

2r 2 

fr(G) < 



2^-2 Q v Q UqV 

Note that uq € V3, since N(v) C V3 for every v £ Vi, as we proved above. Therefore, since f r {G) > 1 



(17) 



co(r) 
„l-5 ) 



it follows by (17) that 





< 


n i-s 




Qv QuqV 


< 


Q vQuqV 


< 



Q vQ uqv 

2r 2 -\- QvQuqv 



2r 2 c (r) < c"'(r) 



,1-6 



,1-5 



(18) 



for an appropriate function c"'(r) of r. 

Note that uq has at least f2(n 1_l5 ) nei ghb ors in V\ U V2 by (fsj and (12), since (f>(n,r) = cu(l) by 



assumption. Thus it follows by |TT| and (14) that ^2 xl£N ( Uo 



3)\{v} dcg(x) 



c" (r)-n 2 



and thus 



ExeN(u„)\iv} fckxl = n ( nl 35 )- Furthermore Q UoV > T,xeN(u )\{v} d5^R b y the definition of Q UQV , and 



thus Qu oV = flln 1 3S ). Moreover Q v Q Uo v = 3S ), since Q v — ^(^)- Therefore it follows by (18) that 



n(n- 3S ) = Q v Q UoV < 

which is a contradiction, since 5 < j by assumption. Therefore there exists no class Q of <7(n)-universal 
amplifiers for any r > ro = 1, where = ^(n 1 ) for some J = \ — e < \. m 

The next corollary follows from Theorem [T] 
Corollary 1 There exists no infinite class Q of undirected graphs which are strong universal suppressors. 

3.2 A class of strong selective amplifiers 

In this section we present the first class Q = {G n : n > 1} of strong selective amplifiers, which we call the 
urchin graphs. Namely, the graph G n has 2n vertices, consisting of a clique with n vertices, an independent 
set of n vertices, and a perfect matching between the clique and the independent set, as it is illustrated in 
Figure [T] For every graph G„, we refer for simplicity to a vertex of the clique of G n as a clique vertex of G n , 
and to a vertex of the independent set of G„ as a nose of G n , respectively. We prove in this section that the 
class Q of urchin graphs are strong selective amplifiers. Namely, we prove that, whenever r > ro = 5, the 
fixation probability of any nose v of any graph G n is f r (v) > 1 — where c(r) is a function that depends 
only on the fitness r of the mutant. 

Let v be a clique vertex (resp. a nose) and u be its adjacent nose (resp. clique vertex). If v is infected and u 
is not infected, then v is called an isolated clique vertex (resp. isolated nose), otherwise v is called a covered 
clique vertex (resp. covered nose). Let k £ {0,1,..., n}, i £ {0, 1, 2, . . . , n — k\, and x € {0,1,2,..., k}. 
Denote by Q\ x the state of G n with exactly i isolated clique vertices, x isolated noses, and k — x covered 



noses. An example of the state Qi x ^ s illustrated in Figure 2(a) Furthermore, for every k, i € {0, 1, . . . , n}, 
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Figure 1: The "urchin" graph G n with 2n vertices. 



we define the state Pf of G n as follows. If i < fc, then Pf is the state with exactly i covered noses and 
k — i isolated noses. If i > k, then Pf is the state with exactly k covered noses and i — k isolated clique 
vertices. Note that Q% — P^ +i and Q§ x — P£_ x , for every k <E {0, 1, ... , n}, i £ {0,1,2, ... ,n — k}, and 
x € {0, 1, 2, ... , k}. Two examples of the state P/°, for the cases where i < k and i > k, are illustrated in 
Figures [2(b)| and 2(c) respectively. 





Pi-. 




(a) (b) (c) 

Figure 2: The state (a) Q\ x and the state Pf, where (b) i < k, and (c) i > k. 

Let k G {1, 2, . . . ,n — 1}. For all appropriate values of % and x, we denote by q\ x (resp. pf) the probability 
that, starting at state Q\ x (resp. P*) we eventually arrive to a state with fc + 1 infected noses before we 
arrive to a state with k — 1 infected noses. 

Lemma 2 Let fc € {1,2, . . . ,n — 1}. Then q^ x > q^_i x _i, for every i £ {1,2, . . . ,n — fc} and every 
are {1,2,..., A;}. 

Proof. Denote by Mi the Markov chain with starting state Q\ x . Similarly, denote by M.2 the Markov chain 
with starting state Qi— i x -i- Note that, initially, both Markov chains M\ and M.2 have the same number 
2(fc — x) + x + i of infected vertices. Moreover, Mi and M2 coincide initially on all their vertices except two. 
In particular, Mi has initially an isolated nose u, which is a covered nose in M2- Furthermore, Mi has 
initially an isolated clique vertex v, which is an uninfected clique vertex in M2- Denote by v! the (unique) 
clique vertex that is adjacent to u in G n . Furthermore denote by v' the (unique) nose that is adjacent to v 
in G n . Note that, initially, v! is uninfected in Mi and infected in M2, while v' is uninfected in both Mi 
and M 2 - 

Note that at every iteration of the processes Mi and M2, one vertex w is activated and then it replaces 
a neighbor w' of it by an offspring of w. Thus, an equivalent way to analyze these processes is to consider 
that, at every iteration, one directed edge between two adjacent vertices is activated (with the appropriate 
probability). In order to prove that qf x > qf_ 1 a ,_ 1 , we simulate the progress of Mi by the random choices 
made at the corresponding steps by M2- In particular, we simulate the processes Mi and M2 until they 
reach states 5i and S2, respectively, such that either Si = S2, or one of Si and S2 is strictly included in 
the other. Furthermore, during the whole simulation of Mi by M2, before we reach such states Si and S2, 
both Mi and M2 have the same number of infected vertices at the corresponding iterations. 
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If a w ^ {u, u', v, v'} is activated for reproduction in M2, then we activate w also in Mi- If w places its 
offspring at a vertex w' ^ {u,u f ,v,v'} in M2, then w places its offspring at the same vertex w' also in Mi- 

If the clique vertex v (resp. u') is activated for reproduction in M2, then we activate u' (resp. v) in Mi- 
In this case, if v (resp. u') places in M2 its offspring at a clique vertex w ^ u' (resp. w ^ v), then v! (resp. v) 
places in Mi its offspring at the same clique vertex w. If v (resp. u') places in M2 its offspring at the clique 
vertex u' (resp. v), then u' (resp. v ) places in Mi its offspring at the clique vertex v (resp. u'); in this case 
we arrive to two identical states in both Mi and M2- Finally, if v (resp. u') places in M2 its offspring at 
its adjacent nose v' (resp. it), then it' (resp. v) places in Mi its offspring at its adjacent nose u (resp. v'); in 
this case we arrive in Mi to a state, in which the infected vertices are a strict subset (resp. superset) of the 
infected vertices inX 2 . 

If a clique vertex w ^ {v,u'} is activated for reproduction in M2, and if w places in M2 its offspring at 
v (resp. it'), then w places in Mi its offspring at v! (resp. v). In this case, if the number of infected vertices 
changes in M2, then we arrive to the same state in both Mi and M2- 

Finally, if the nose w — v' (resp. w — u) is activated for reproduction in M2, then we activate the same 
nose also in Mi- In this case we arrive in Mi to a state, in which the infected vertices are a strict subset 
(resp. superset) of the infected vertices in M2- 

Note now that q\ x > q\_ x x _ 1 if and only if, in the above simulation of Mi by M2, the probability 
that we arrive to strictly more infected vertices in Mi than M2 is greater or equal to the probability that 
we arrive to strictly less infected vertices in Mi than M2- Furthermore, note that we arrive in Mi to a 
state with strictly more or strictly less infected vertices than in Ai 2 only when one of the edges uv! or vv' 
is activated (in some direction) in the process M2- In particular, whenever this event occurs, Mi receives 
strictly more infected vertices than M2, if either u places its offspring at u' in M2, or if u' places its offspring 
at u in M2- Similarly, Mi receives strictly less infected vertices than M2, if either v places its offspring at 
v' in M2, or if v' places its offspring at v in Mi. The ratio of these probabilities is 

r ■ - + r ■ 1 

? = r > 1 

- + 1 

n 

and thus q k i x > q%_ 1)X _ v ■ 

Corollary 2 Let k € {1, 2, . . . , n - 1}, i € {0, 1, . . . , n - k}, and x € {0, 1, ... , k}. Then q\ x > p\ +i _ x . 

Proof. Suppose first that i > x. Then Lemma [2] implies that qf x > q^i^i > ... > ?|L X! o = Pk+i-x- 
Suppose now that i < x. Then Lemma [2] implies that q\ x > <^_x,x-x > ■ ■ ■ > 1o,x-i = Pk+i-x- ■ 

Note by Corollary [2] that, in order to compute a lower bound for the fixation probability f r {v) of a nose v 
of the graph G„, we can assume that, whenever we have k infected noses and i infected clique vertices, we 
are at state P i K That is, in the Markov chain of the generalized Moran process, we replace any transition 
to a state Q\ x with a transition to state P^ +i _ x - Denote this relaxed Markov chain by M; we will compute 
a lower bound of the fixation probability of state Pq in the Markov chain M (cf. Theorem [3]). 

In order to analyze M, we decompose it first into the n — 1 smaller Markov chains Mi, M2, ■ ■ ■ , M n -i, 
as follows. For every k 6 {1, 2, . . . , n— 1}, the Markov chain Mk captures all transitions of M between states 
with k infected noses. The state graph of Mk is illustrated in Figure|3| where we denote by F k _i (resp. F k+ i) 
an arbitrary state with k — 1 (resp. k + 1) infected noses. Moreover, we consider Fk-i and -Ffc+i as absorbing 
states of Mk- Since we want to compute a lower bound of the fixation probability, whenever we arrive at 
state Fk+i (resp. at state Fk-i), we assume that we have the smallest number of infected clique vertices 
with k + 1 (resp. with k — 1) infected noses. That is, whenever Mk reaches state -Ffc+i, we assume that M 
has reached state P^i (and thus we move to the Markov chain Mk+i)- Similarly, whenever Mk reaches 
state Fk-i, we assume that M has reached state Pg _1 (and thus we move to the Markov chain Mk-i)- 

3.2.1 A decomposition of Mk into two Markov chains 

In order to analyze the Markov chain Mk, where k £ {1, 2, . . . ,n — 1}, we decompose it into two smaller 
Markov chains {M\,M\}, as they are shown in Figure El 
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F k -i 



Fk+i 



F k +i 



I I 



Figure 3: The state graph of the relaxed Markov chain A4k, where k £ {1, 2, . 



!}• 



J I) 



J 2 






Fk-i 



JDk 
r k+l 



Figure 4: The two Markov chains Ai\ and Mf,, where k £ {1, 2, . . . , n — 1}. 
In M\, we consider the state P^ +1 absorbing. For every i £ {0, 1, . . . , k} denote by h\ the probability 



that, starting at state in M.\, we eventually reach state P^ +1 before we reach state -Ffc-i, cf. Figure 4(a) 



In this Markov chain Ai\, every transition probability between two states is equal to the corresponding 
transition probabilities in Aik- 

In A4fr, we denote by sf, where i £ {k,k + l,...,n}, the probability that starting at state Pf we 



eventually reach state Ff. + i before we reach state Ff._i, cf. Figure 4(b) In this Markov chain .MS, the 



transition probability from state P£ to state P£ +1 (resp. to state Fk-i) is equal to h\ (resp. 1 — h\), while 
all other transition probabilities between two states in A4f, are the same as the corresponding transition 
probabilities in Aik- 



Lemma 3 In the Markov chain M.\, for any r > 1, 

h k k >l- 2 =1-0(1) 
n(r — 1) + 1 n 

Proof. For i = 0, the value of h\ in the Markov chain M\ is 



(19) 



rk ■ h\ + | • 



rk 



rn + 1 



■hi 



(20) 
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and thus 



h\ - h\ 



11 rn rn 



(21) 



Furthermore, for every i £ {1. . . . , k}, where 1 < k < n — 1, the value of /i^ in A'iJ, can be computed as 
follows. 

h k = a k i -h k i+l +^-h k _ t + 1 k -Q (22) 

where 



i ((fc— i)) 



O ; 



# = 



i(n— i) 
n 

E* 

k—i 
n 

E? 



(23) 



and E, fc = r ((fc-»)n+»(»-^)) 



i(n-i) k—i 



Therefore (22 1 implies that 



4(^-^-1) + 



(24) 



Furthermore, (23) implies that 



01 



1 



< 



1 



r (k — i)n + i(n — i) r 
1 k-i 1 



r (k — i)n + i(n — i) r n \ H"— 



< 



(25) 
(26) 



Note that the inequality < ~ in (26) holds also for i — k, since 7^ = 0. Therefore it follows by (24), (25), 
and ((26} that 



n i+i 



h k <-(h k -hl 1 ) + -h k 
r rn 



Thus, since h k ^ — 1 by definition, it follows by |27| for i — k that 



(27) 



1-^ < -(h k -h k _ 1 ) + -ii 

r rn 



1 



- K-2) + 



1 



1 1 



l fc-i 



< 

< ... 

< \(h k h k ) + -(h k + -h k _ x + ... + -]^h k ) 



Therefore, since h\ < h\ < . . . < h\ and h\ — /ig < ~ by (21), it follows that 

1 



,.A- + 1 



- + -^(1 + ^ 

n rn r 



1 , 1 

< 



„k-l 



'hi 1 



n rn 1 — 



1 1 fc < 1 1 fc 

„fc+l n n ( r _l) fe -(r-l)?i n(r-l) k 



Therefore 



n(r — 1) — 1 
ri(r - 1) + 1 



11 



i.e. 



h k k >l 



n(r - 1) + 1 



Lemma 4 In the Markov chain for any r > 1, 



fc + 2 



n(r — 1) 

Proof. Let 1 < i < k. Recall by (27) in the proof of Lemma [3] that 



< 




-hU) 


1 i k 

+ —h- 
rn 


< 








< 




- K) 4 


rn r 


< 






1 , 1 
--(1 + - + .. 
rn r 


< 




- h k ) 4 


1 

n(r — 1) 



— h k ) 



Therefore it follows by (21 1 that 



n i+l 



hf < 



1 



r' Ti n n 



(r-1) 



Summing up (30) for every i = 1, 2, . . . , k, it follows that 

1 - h\ < 



1 - 1 



1 



< 



n r^ 
1 1 



.fc+l 



n(r — 1) 
fc + l 



< 



since h\ +1 = 1, and thus 



n r[r — 1) n(r — 1) n{r — 1) 
fc + l 



hi > 1- 



n(r — 1) 



Therefore it follows now by ( 20 ) that 



W > 



"(1 



fc+l n(r - 1) fc + l 

' ( ~~ n(r- l) j 

fc + 2 



rn + 1 71(7' — 1) n(r — 1) + 1 
n(r — 1) — fc — 1 
n(r — 1) + 1 ~~ 



fc + 2 



> 1 



n(r — 1) + 1 n(r — 1) 



We now provide some lower bounds for the Markov chain M\. 
Lemma 5 In the Markov chain M\, for any r > 5, 

4 > 1 



(r — 5)(r-l) (n-fc) 2 
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Proof. For i= k, the value of s k in the Markov chain M\ is 



sl = ht-st +1 + {l-hi)-Q = h k k -s k 



Therefore Lemma [3] implies that 

s k k >(l 

and thus 



n(r - 1) + 1 



)4+i > (i 



k °k+l 

2 

n(r — 1) 
2 



Sfc+1 ~~ Sfc " n(r-l) Sfc+1 " n(r - 1) 



Furthermore, for every i G {fc + 1, . . . , n}, the value of s k in M.\ can be computed as follows. 



where 



a k = 



k -s k +1 +p k -s k _i+J--l 



i(n—i) 



(33) 

(34) 
(35) 

(36) 



I 

U — JA -L ^( n ~0 (t-fc)n+i(n-i) 
>fc V n n 



Pi = 



7* = 



(37) 



, i—k 



and E k = r 1 -^- + (t ~ fc) "+ t( "~' ) + r*=*. Therefore g implies that 



Furthermore, (37) implies that 



CV; 



(i — fc)n 
i(?i — i) 



i — k i — k 1 
i(n — i) ~ i n 



(38) 

(39) 
(40) 



We now prove that §s whenever i < =±*. Suppose first that fc < f . Then i < ^ < ^±^, i.e. i < ^. 
Thus < and thus (39) implies that ^ < H 1+4 ) = f . Suppose now that n > k > § . Then also i > § , 
since £ > k + 1, and thus - < 2. Furthermore i-k<n-i, since z < Therefore = • ? 

and thus (39) implies that < £ • (1 + 2) = 2. Summarizing, for every fc £ {1,2 



i6{fc + l 

Therefore it follows by @, and (|4l]) that 



# < 5 
a k ~ r 



r 



i — k 
in 



^7 — = — • a < 2, 
— 1} and every 

(41) 
(42) 
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Thus, in particular 



sti < -(*f_i-*?_ 2 ) 



(43) 



Now (42) and (431 imply that 



h+i - st < (-r fc (^ +1 - Sk) 



i — k 



(i -4) 



(44) 



Note that s* = S £ + ~ s t) + ( s fc+2 " 4+i) + ■ ■ ■ + (4 ~ 4-i)- Thus @ implies that 

< 4 + (s k k+1 -s k k )-(l+*+... + (*) i - k - 1 ) 



< 



4 + ( s *-+i 4) 
4 + ( s A-+i — 4) 



i 



i - 



r — 5 



(45) 



Therefore ( 44 ) and ( 45 ) imply that 



4+i-4<(-Y- k (4+i-^ 



—(i-4-(4 + i-st)-^ T ) 

m r — 5 



(46) 



Note that ( 46 ) holds also for i — k and that in this case it becomes an equality. Summing up ( 46 ) for every 

2 



i G {k, . . . it follows that 



5. n+fc 



fe+i — s fc) 



E 



i — k, 



< — - -S) - a - -S - -4)--^) E — 



r — 5 



r — 5 * — ' zn 

i=fc+i 



Note now that for any positive numbers x,y,z,w > 0, it holds that - + ^ > 
" {*H 1 }• 



r — 5 * — ' in 

i=fe+l 



(47) 

Therefore, for every 



Thus 



i-k m^-i + k + l)-k 
in (2±* - j + _|_ i) n 



> 



- Jfe + l 
(^ + fc + l)n 

n — k + 2 n — k 

> 



n(n + 3fc + 2) n(n + 3fc) 
n — k (n — k) 2 



i=fc+l 



n(n + 3fc) 2n(n + 3fc) 



E 

«=fe+i 



i — k (n . — k) 2 



in An{n + 3fc) 



(48) 
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It follows now by ( 47 ) and ( 48 1 that 



o < si +k 



Sft 



< 



r — 5 



/ k k\ ^) /-i k l k k\ ^ \ 

[Sk+1 " Sk) " 4n(n + 3fc) (1 " Sfe " (Sfc+1 " Sfc) ■ 7^5 ' 



r — 5 



(n — fc) 2 (n — fc) 2 



/ ft 

l s ft+i s ft 



4n(n + 3fc) 4n(n + 3fc) 
fe 4n(n + 3fc) + (n - fc) 2 (n - fc) 2 



(1 -4) 

11-^) 



Therefore 



and thus 



r-5 v ~ K+1 ~ K/ 4n(n + 3fc) 4n(n + 3fc) 

(n - fc) 2 (l - s k k ) < -\(s k k+1 s k k )(An(n + 3k) + (n - fc) 2 ) 
, ju i.. 4n(n + 3fc) . 



r — 5 



(n-fc) 2 



Now ( 49 ) and ( 35 1 imply that 



4 > i- 
= i - 



?■ 



2 | 4n(n + 3fc) ) > 1 



r — 5 n(r — 1) v ' (n — fc) 2 
16r n + 3fc 



2 2 4n(n + 3fc) 



> 1 - 



r — 5 n(r — 1) (n — fc) 2 
64r n 



(r-5)(r-l) (n-k) 2 ~ (r - 5)(r - 1) (n - fc) 



(49) 



The next two corollaries follow now from Lemma [5] by direct substitution. 
Corollary 3 In the Markov chain M\, for any r > 5 and any fc < f , 

i- 64r 4 ^.1, 

s fc > 1 - 7 ^7 7T ' " = 1 - O (-) 

fe ~ (r - 5)(r - 1) n n' 

Corollary 4 In the Markov chain M\, for any r > 5 and any k < n — \Jn logn, 



4 > 1 



^ L_ = i_ 0( JL, 

(r — 5)(r — 1) logn logn 



We now present an auxiliary lemma that provides a trivial lower bound for the probability s k , for any 
fc < n — 1. 

Lemma 6 In i/ie Markov chain A4 2 , for any r > 5 and any fc < n — 1, 



1 



Proof. Let 1 < k < n— 1. Recall by (36) and (37) in the proof of Lemma [5] that for i = fc + 1 



■ 4+2 + Pk+i ■ 4 + ik+i ■ 1 



'fc+i 



Therefore, since s^ +2 > s k+1 , it follows that 



ft Jtr4+7tri 



$+1 + 7^+1 



(50) 
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In particular, it follows by (50) and (371 for i = k + 1 that 

(n + (k + l)(n - k - 1)) • s£ + r ■ 1 



s fc > 



n + (fc + l)(n - A; - 1) + r 
Furthermore recall by ( 34 ) in the proof of Lemma [5] that 

n(r-l) u 2 



n(r — 1) — 2 n(r — lj — 2 



(51) 



(52) 



Thus ([51) and (|52j) imply that 
(1 



fc (n+(fc + l)(n-fc-l))-«£+r-l 



n(r - 1) - 2 



'*fc> 



n+(k + l)(n - fc - 1) + : 



and thus 



(n+(fc + l)(n- A; - 1)) 



(n + (k + l)(n - A; - 1) + r)(l + 
2 



n(r - 1) - 2 



> (n + (A; + l)(n - fc - l))s| + r ^ 



n(r-l)-2~ R ' ' n(r-l)-2 / " fe 
(2(n+(fc + l)(n-fc-l))+rn(r-l))s£ > r(n(r - 1) - 2) 



Therefore 



s k — 



r(n(r- 1) - 2) 



2(n + (fc + l)(n - A; — 1)) + rn(r - 1) 
Note now that (k + l)(n — A; — 1) < n 2 , and thus the last inequality implies that 



fe ^ r(n(r-l)-2) > 



4> 



(r - 1) - - 



2(n + n 2 ) + rn(r - 1) ~~ 2(n + 1) + r(r - 1) 
Therefore, since r > 5 and r < n by assumption, it follows that 



4> 



20-2 



9 1 

> 



2(n + l) + 20 n + 11 n 



3.2.2 Urchin graphs are strong selective amplifiers 



In this section we conclude our analysis by combining the results of Section 3.2.1 on the two Markov chains 
M.i and M. 2 . The Markov chain M. is illustrated in Figure [5J where the transition from state Pq to the states 
P£, Pq _1 is done through the Markov chain M.\, and the transition from state P£ to the states P^+ii Pq -1 
is done through the Markov chain M.i 1 respectively. 






Figure 5: The Markov chain M., using the Markov chains M\ and M.\, where k e {1, 2, . 
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In Figure [HJ the transition probability from state P£ to state P^+l (rcsp. Pq 1 ) is (resp. 1 — s\). 
Recall that s\ is the probability that, starting at P^ in M. 2 (and thus also in Ai), we reach state F^+i before 
we reach Fk-i- Furthermore, the transition probability from state Pq to state Pj? is equal to the probability 
that, starting at Pq in Ai%, we reach P* before we reach Ffc_i. Note that this probability is larger than Iiq. 
Therefore, in order to compute a lower bound of the fixation probability of a nose in G n , we can assume 
that in Ai the transition probability from state Pq to P^ (resp. Pq _1 ) is (resp. 1 — h^), as it is shown in 
Figure [5] 

Note that for every k € {2, . . . , n — 1} the infected vertices of state Pq is a strict subset of the infected 
vertices of state P£. Therefore, in order to compute a lower bound of the fixation probability of state Pq 1 



in Ai, we can relax Ai by changing every transition from state P k Z\ 



to state P^ to a transition from state 

" After 



1}. This relaxed Markov chain Ai' is illustrate in Figure 6(a) 
, n 



Pfcli to state P fc , where k € {2, . . . 

eliminating the states P^ in Ai' , where k g {1, 2, . . . , n — 1}, we obtain the equivalent birth-death process 
B n that is illustrated in Figure 6(b) Denote by p\ the fixation probability of state Pq 1 in B n , i.e. p\ is the 
probability that, starting at state Pq 1 in B n , we eventually arrive to state P™. 





jk-l 1 fy) pk 
n r n 




(a) 



1 — ft 1 * 1 




P l hls\ PS 



i - K4 





Figure 6: (a) The relaxed Markov chain Ai' and (b) the birth-death process B n that is obtained from Ai' 
after eliminating the states P^ in Ai' , where fee {1, 2, . . . , n — 1}. 



Theorem 2 Por any r > 5 and for sufficiently large n, the fixation probability p\ of state Pq in B n is 
Pi > 1 — ^p- , for some appropriate function c(r) of r. 

l h k s k 

Proof. Denote by Afc the forward bias of B n at state Pq , i.e. = k k k is the ratio of the forward over 
the backward transition probability at state Pq . Then the fixation probability p\ of state Pq 1 in B n is 

Pl = 1 + ^ + ^ + ^ + ... + (53) 

Ai A1A2 A1A2A3 Ai A2 A3 . . . A^ — 1 



Note now by Lemma [3] and Corollary [3] that for every k 



< 



uk k (\ - k + 2 \(-\ - 256r . 1\ 

a "oSfc ^ ^ n(r-l)^ 1 (r-5)(r--l) n > 
At- = n— r > 



,fc„fc — -1 _ f-i _ fc+2 \fi _ 256r 1\ 
( 0*fc 1 I 1 nfr-l)^ 1 (r-5)(r-l) ' 



1_ 256r 1_ fc+2 n(r _ 1} _ 256r _ (fc + 2 ) 

(r- 5)(r-l) n n(r-l) _ V / (r-5) V / ,_ ... 

— 256 r 1 I fc+2 256r i fj. i o\ \ D ' 

(r-5)(r- 1) ' n + n(r-l) 0"-5) " r ^ " r ^ 
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Therefore, since k < § and 



and thus 



256r 
(r-5) 



< logn < | — 2 for sufficiently large n, it follows by (54) that 



logn + (fc + 2) 21ogn + fc 



1^2 log n + k 



\ k n(r-2) 

for every k < ^. Furthermore, note by Lemma [4] and Corollary [4] that, whenever j < k < n — \Jn logn, 

(1 _ fc+2 Wl _ 64r 1 \ 

^ n(r-\)>y L (r-5)(r-l) log n ^ 



(55) 



A, 



0*k 



> 



> 



1 _ Uk k — 1/1 fc+2 \fi 64r _ 1 

1 "0*fc 1 I 1 n(r-l) A 1 (r-5)(r-l) logrJ 

1- . g£ n ■ ^ £±2 (r-l)logn- ^logn 

(r— 5)(r— 1) log n n(r— 1) V y o ( r_ 5) n to 



64r 



. 64r , fc+2l 

(r-5)(r-l) logn ^ n(r-l) (r-5) ' n lu & " 



1 , fc+2 



Therefore, since k + 2 < n and 



(r-5) < 7=^ f° r sufficiently large n, it follows that 
(r-l)logn-^-log 

log n 
r-2 



^ ^logn 



r-3 



and thus 



1 1 

— < 

A fc r-3 



whenever ^ < k < n — \Jn log n. Moreover, note by Lemma [2] and Lemma [6] that, whenever n — \Jn log 
k < n — 1, 



(56) 
n < 



k c k 



h!ks 



0°k 



> 



(I _ \ 1 

v n(r— 1) ' n 



1 - 



i(r— 1) ' 

1 _ n+1 i 2 



> 



n n 



>(r-l) 



> 



1- 1 + 2 ? ? +1 n n-l + ^-r 
n ' n J (r— 1) r— 1 



Note now that 1 — -^j > | since r > 5 by assumption, and thus the latter inequality implies that A^. > i^-, 
i.e. 



A ft 



< 2n 



(57) 



whenever n — \Jn log n < k < n — 1. 



Since r > 5 by assumption, note now by (55) that ^- < ^('"-si) wnenever ft < logn, and that 



X^T < n (r-2) ^ n whenever logn < fc < Therefore, for every k £ {2,3, ... , logn}, 



1 



A1A2A3 ... Afc 
Furthermore, for every k € {log n + 1, . . . , §}, 

1 



< 



Ai A2A3 ... Aj, 



< 



< 



log n 



log n 



log n 



< 



log n 



n 

i— log n+1 



< 



logn 



(58) 

(59) 
(60) 
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Therefore, for every k £ {2, 3, . . . , 



E 

fc=2 



A1A2A3 . . . Afc 



< logn 



2 , \ 3 
logn\ n /logn 



log d n log d n _ 3 log 13 n 1 
n 2 2n 2 2n 2 n 



(61) 



for sufficiently large n. Note furthermore by (56) that j- < 1 whenever ^ < k < n — -^/nlogn, since r > 5 
by assumption. Therefore, for every fce{^ + l,...,n — ^/nlogn}, 



< 



A1A2A3...AJ; A1A2A3 . . . A2„/3 



and thus it follows by ( 60 ) that 



n—^/n log n 

E 



1 



A1A2A3 . . . Af, 



< 



n log n 



logn\ log n 1 

< 



2n 2 n 



for sufficiently large n. Let now n — \J n log n < k < n — 1. Then it follows by (57) and (59) that 

log n 



1 



A1A2A3 ... Afc 



< 



< 



log n 



log n 



i— log 71+1 



n -,-^) k 



k — n-\-y/r\ 



log n 



(o) 



i/2-logn _ (2^)^" log ™ 



2 los ™ (logn) log ™ ■ 2^™ logn • n^ nlogn 



ilogn 



2 n/2 



< 



2log n log log n _ 2 V n l°g n . 2 ^°S n V n l°g n 
2W2 

2(log n log log n-\-y/n log n-flog ny/n log 
2^2 



However 



log n log log n + \J n log n + log ny^n logn < 



for sufficiently large n, and thus (63) implies that 



1 



< 



1 



A; A, A,... A;. 2«/4 

for every k G {n — V n log n + 1, . . . , n — 1}. Therefore 

n-i 



E 



fc=n— Vn log n+l 



1 n 1 

AiA 2 A 3 ... A fc < 2^ < n 



for sufficiently large n. Thus, summing up (|61j), (|62), and ( |64[ ), it follows that 

1 



n-l 

E 

fc=2 



3 

AiA 2 A 3 ...Afe n 



(62) 



(63) 



(64) 



(65) 
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For k = 1, (54) implies that 



A fc > 



n(r — 1) 



(r-5) 



256r 
(r-5) 



> n- 



(r-2) 

256r j_ <: 



(r-5) 



(r-2)(r-5) 
259r - 15 



and thus 



< 



259r - 15 
(r - 2)(r- 5) 



Summarizing, it follows by (53), (65), 



Pi = 



1 

and (|66j that 
1 



(66) 



259r-15 



> 1 - 



c(r) 
n 



(r-2)(r-5) 

where c(r) is an appropriate function of r. This completes the proof of the theorem. ■ 

We are now ready to provide our main result in this section. 

Theorem 3 The class Q = {G n : n > 1} of urchin graphs is a class of strong selective amplifiers. 

Proof. Consider the urchin graph G„, where n > 1. Let v be a nose in G n . Then the fixation probability 
f r (v) of v in the generalized Moran process is greater than or equal to the fixation probability of state Pq 1 in 
the Markov chain M. of Figurej5] (cf. Corollary [2] and the discussion after it in Section 3.2 ). Furthermore, the 
fixation probability of state Po in the Markov chain M is g reater than or equal to the fixation probability 
Pi of state Pq in the birth-death process B n in Figure 6(b) Therefore, since p± > 1 — for any r > 5 by 

Theorem [2J it follows that f r {v) > 1 — for r > ro = 5 and sufficiently large n, where c(r) is a function 
that depends only on r. Finally, since there exist exactly ^ noses in G„, it follows by Definition [2| that 



the class Q of urchin graphs is a class of ( ^ , n)-selective amplifiers, and thus Q is a class of strong selective 
amplifiers. ■ 



4 Suppressor bounds 

In this section we prove the first non-trivial lower bound for the fixation probability of an arbitrary undirected 



graph, namely the Thermal Theorem (Section 4.1), which generalizes the analysis of the fixation probability 
of regular graphs 16 . Furthermore we present for every function (/>(n), where 4>{ n ) = an d 4>{ n ) ^ V^i 
a class of (^S+T' 07^) )-selective suppressors in Section 4.2 



4.1 The Thermal Theorem 

Consider a graph G — (V,E) and a fitness value r > 1. Denote by M. r {G) the generalized Moran process 
on G with fitness r. Then, for every subset S £ {0, V} of its vertices, the fixation probability f r (S) of S 
in M r {G) is given by ([T]), where / r (0) = and f r (V) = 1. That is, the fixation probabilities f r (S), where 
S ^ {0, V}, are the solution of the linear system ([T]) with boundary conditions / r (0) = and f r (V) = 1. 

Suppose that at some iteration of the generalized Moran process the set S of vertices are infected and that 
the edge xy € E (where x € S and y S) is activated, i.e. either x infects y or y disinfects x. Then (JlJ implies 
that the probability that x infects y is higher if is large; similarly, the probability that y disinfects x 



is higher if is large. Therefore, in a fashion similar to 16 , we call for every vertex v £ V the quantity 
d e SV the temperature of v. a "hot" vertex (i.e. with high temperature) affects more often its neighbors than 
a "cold" vertex (i.e. with low temperature). 

It follows now by ([T]) that for every set S ^ {0,^} there exists at least one pair x(S),y(S) of vertices 
with x(S) € S, y(S) S, and x(S)y(S) £ E such that 

dcga:(S) degy(S) 
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Thus, solving the linear system that is obtained from (67 1 by replacing inequalities with equalities, we obtain 



a lower bound for the fixation probabilities f r (S), where S {0, V}. In the next definition we introduce a 
weighted generalization of this linear system, which is a crucial tool for our analysis in obtaining the Thermal 
Theorem. 

Definition 3 (the linear system L ) Let G = (V,E) be an undirected graph and r > 1. Let every vertex 
v 6 V have weight (temperature) d v > 0. The linear system Lq on the variables p r (S), where S C V , is 
given by the following equations whenever S £ {0, V}: 

rd x{S )Pr(S + y(S)) + d y{s)Pr (S - x{S)) 

Pr(S) = — (68) 

ra x{S) + a y{S) 

with boundary conditions p r (0) = and p r (V) = 1. 

With a slight abuse of notation, whenever S — {u\, v,2, . . . , itfe}, we denote p r (ui, U2, . . . , life) = p r {S). 

Observation 1 The linear system Lq in Definition [5] corresponds naturally to the Markov chain A4q with 
one state for every subset S C V, where the states and V are absorbing, and every non-absorbing state S 
has exactly two transitions to the states S + y(S) and S — x(S) with transition probabilities qs = rd ( 
and 1 — qg, respectively. 

Observation 2 Let G = (V,E) be a graph and r > 1. For every vertex x £ V let d x = de ^. a , be the 
temperature of x. Then f r (S) > p r (S) for every S C V, where the values Pr (S) are the solution of the linear 
system Lq. 



Before we provide the Thermal Theorem (Theorem 



lemma which generalizes the Isothermal Theorem of 16 
number of neighbors for every vertex. 



4j, we first prove an auxiliary result in the next 



for regular graphs, i.e. for graphs with the same 



Lemma 7 Let G — (V, E) be a graph with n vertices, r > 1, and d u be the same for all vertices u 6 V. 
Then for every vertex u £ V , 

1-i 1 

Pr{u) = ~ _ [_ > 1 - 



r 



Proof. Since d u is the same for all vertices u G V, it follows by (68 1 that for every set S ^ {0, V}, the 
forward probability is qs = and the backward probability is 1 — qs = ^rpj- Therefore, by symmetry, 
Pr(S) = Pr(S') whenever \S\ = \S'\. For every < k < n denote by P k = Pr{S), where \S\ = k. Note that 



Po = and p n = 1. Then it follows by (68) that, whenever 1 < k < n — 1 



Pk+i ~Pk = -(Pfc -Pk-i) = ••• = — fe (Pi - Po) 
Therefore, summing up these equations for every 1 < k < n — 1 it follows that 



w 1 1 1 N 

p n -pi = (pi - po){- + - 2 + ■ • • + -^r) 



and thus, since p Q = and p n = 1, 



1 _ 1-^ 1 

Pl ~ i + i + ^ + ... + ^r ~ r^fc - 1 ~ r 



We are now ready to provide our main result in this section which provides the first non-trivial lower 
bound for the fixation probability on arbitrary graphs, parameterized by the maximum ratio between two 
different temperatures in the graph. 
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Theorem 4 (Thermal Theorem) Let G = (V,E) be a connected undirected graph and r > 1. Then 
fr(v) > ■ r dc S 5 for every v € V. 

T+ dcg min 

Proof. Let G have n vertices, i.e. |V| = n. Our proof is based on the linear system Lq of Definition [3] 
Namely, we consider the linear system L with weight d v = for every vertex v £ V. Note that 

rfmin = — and d max = — . Recall that Aio is the Markov chain that can be defined from the linear 
system L$ (cf. Observation [jj, and that every state 5 ^ {0,1^} of A^o has exactly two transitions, namely 
to states S + y(S) and S - x(S). 

We now define the Markov chain A4q from Ai as follows. Consider an arbitrary state S C V such that 
l < \S\ < n — 2. Denote x(S) — u and y(S) = v (note that v ^ S). Furthermore denote x(S + v) = x and 
y(S + v) = j/o- Then perform the following changes to the Markov chain A4o : 

Step A. add a new dummy state Xs to Aio, 

Step B. replace the transition from S to S + v by a transition from S to (with the same transition 
probability qs), 

Step C. add to state the transitions to states S + v + yo and S + v — xq, with transition probabilities 
qs+v and 1 — qs+v , respectively. 

An example of the application of the above Steps A, B, C is illustrated in Figure [7] Denote by A4q 
the Markov chain obtained after applying these steps to M. for every state S C V with 1 < \S\ < 11 — 2. 
Furthermore denote by L'q the linear set that corresponds to Mq- Note that Lq has the additional variables 
{p r (Xs) ■ 1 < \S\ < n — 2} that do not exist in Lq. Moreover, for every state Xs of Mq, note by the 
construction of A^g that p r (Xs) = p r (S + y(S)) in the solution of Lq. 




Figure 7: Parts of (a) the Markov chain A4q and of (b) the Markov chain A4q (after the execution of 
Steps A, B, and C. 

Consider an arbitrary numbering vq, v\, . . . , v n —i of the n vertices of G. For every i £ {1, 2, . . . , n — 1}, 
we construct iteratively the linear system L* (and the corresponding Markov chain M.*) from L*_ l (and 
from the corresponding Markov chain M.*_-j) as follows. Consider a state S £ {0, V}, where y(S) — Vi (note 
that in this case Vi (£ S). Denote x(S) — u, x(S + Vi) — xq, and y(S + Vi) = yo. Then perform the following 
changes to the Markov chain M.*_ x : 

Step 1. replace the transition from Xs to S + Vi + yo by a transition from Xs to S + yo (with the same 
transition probability qs+vi), 

Step 2. replace the transition from Xs to S + Vi — xo by a transition from Xs to S — xo (with the same 
transition probability 1 — qs+vi ) ■ 

Denote by Ai* the Markov chain obtained by iteratively applying Steps 1 and 2 to for every state 

S {0, V} with y(S) = u,. Note that for every state Si ^ {0, V} and for every state S% 7^ V, where Vi G S2 
and Vi Si, there exists no transition path in Ai* from Si to 62- Furthermore denote by L* the linear 
system that corresponds to Ai* (cf. Observation [I]) . Note that, in the above Step 2, if xo = ^ 5 then 
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V l-qs 



qs+v. 




gs+ 



V l-qs 



• S - u • S — u 

(a) (b) 

Figure 8: Parts of (a) the Markov chain and of (b) the Markov chain Ai*, where xo = Vi 

S X s S + v, + y S X s S + y 



VI- q s 



S — u 



qs+v, 
i - qs+v, 

• S + Vi - Xq 

(a) 



'Is 



VI -q s 



S — u 



qs+v, 
i - qs+v, 



s - x 



(b) 



Figure 9: Parts of (a) the Markov chain Ai*_ 1 and of (b) the Markov chain Ai*, where Xo ^ Vi 



S — xo — S, and thus we have in this case a transition from Xs to S in M*. The changes made to M*_ 1 by 
Steps 1 and 2 are illustrated in Figures [8] and [9] for the cases where x — v t and x ^ respectively. 

For every state S (resp. Xs) of the Markov chain M*, and for any i 6 {0, 1, . . . , n— 1}, we denote in the 
following by p l r {S) (resp. (Xs)) the value of p r (S) (resp. p r (Xs)) in the solution of the linear system L*. 
Note that for every state S {0, V} with y(S) = Vi, the transitions of state Xs in A4*_ 1 are to the states 
S + Vi + Do and S + Vi — xo, while the transitions of state Xs in Ai* are to the states S + yo C S + v t + y a 
and S — x <Z S + Vi — xq, respectively. Therefore p l r (Xs) < p t r ~ 1 (Xs), and thus p l r {S) < p t r ~ 1 (S). Similarly 
P l r {T) < p*~ 1 (T), for every state T CV with Vi ^ T. Thus in particular p l r (vj) < p* _1 (w J ) for every j ^ i. 

In order to continue our analysis, we distinguish now the cases where xq = Vi and xq Vi. 

Case 1: x a — v i: cf. Figure [Hj In this case we have in the linear system L* that 

pl(S) = qspl.(Xs) + (l-qs)p i r (S~u) 
p*.(X s ) = qs+^pKS + y ) + (1 - qs+vMiS) 



and thus 



where 



pUS) 



<!SqS+v t P l r(S + yo) + (1 - qs)p l r (S - u) 

qsqs+vi + (i - qs) 



qs 
qs+vi 



rdu 
+ 
rd v 



rd,. 



Therefore the forward probability of state S in Ai* is (after eliminating the state Xs) equal to 

qsqs+v, _ r 2 d u 

qsqs+vi + (1 - qs) r 2 



(69) 
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Case 2: Xq ^ Vi, cf. Figure [9j In this case we have in the linear system L* that 

pUS) = qspUXs) + (1 - qsH(S - u) 

P l r( X s) = qS+v z P l r ( S + Vo) + i 1 - 1S+v,)p l r(S ~ X ) 

and thus 

Pr(S) = QSQS+vtPr{S + Vo) + 9s(l - QS+vJpKS - Xq) + (1 - qs)p l r (S - u) 

where 



qs 

qs+vi 



rd u 



rd u + d Vi 

rd Xo 
rd, 4- d 



Therefore the forward probability of state S in Ai* is (after eliminating the state Xs) equal to 

Qsqs+vi = , — -, -, — — i— (70) 

TOj U ~r C*Ui T&xn l 



It follows now by Cases 1 and 2 (cf. (|69|) and (70)) that the forward probability of state S in Ai* (after 



eliminating the state Xs) is a monotone decreasing function of d Vi . Therefore, for every state S' C V with 
Vi ^ S' , the value is also a monotone decreasing function of d Vi . Thus, in particular, also the value 

p l r (vj), where j ^ i, is a monotone decreasing function of We now increase the value of d Vi to ci max in 
L*. Thus for every j ^ i, the value p l r (vj) decreases after this change. 

Recall that f r (vo) > Pr( v o) by Observation [2J Therefore, since also p°(i>o) > ■ • ■ > p"~ 2 ( w o) > p"~ 1 ( u o), 
it follows that f r (vo) > P™ 1 (i'o), i-e. (vo) is a lower bound for the fixation probability f r (vo) in the 
Markov chain A^ r (G). Furthermore d Vl = d V2 = . . . = d Vn l = d max in the linear system L* n _ 1 . Consider 
now the state 5 = {«o} m the Markov chain M-n-ii an d ^ y{S) — Vi , where 1 < io < n. — 1. Note that 
= fo- Then the value p" _1 (i>o) equals 



r4 p™ HwOjVio) + d v t0 Pr *( 

rd VQ n-l 



^K,^) (71) 



. ^ B «-v- i 



cf. Definition jSJ Recall that d Vo — dc ^ Vo and d max = — by definition. Thus, since f r (v ) > p" 1 («o) as 



we proved above, ( 71 ) implies that 



r 



f r (v ) > P n r -\v a ) > tf-i( Vio ) (72) 

Now, similarly to the above transformations of the linear system L*_ 1 to L*, where 1 < i < n—1, we construct 
the linear system L* (and the corresponding Markov chain M*„) from L* l _ 1 (and from the corresponding 
Markov chain M*_ x ), by applying iteratively the above Steps 1 and 2 to the states S C V, where y(S) = v a 
(instead of y(S) = Vi above). Furthermore we increase the value of d VQ to d max in the resulting linear system 
L*. Then, similarly to the construction of L*, where 1 < i < n — 1, it follows that p"(Uj) < p™ _1 (vj) for 
every j ^ 0. Thus, in particular, p"(Uj ) < p™ _1 (ui ). Furthermore d Vo = d Vl = . . . = d v * = rf max in L 
and thus p™(i>i ) > 1 — ^ by Lemma Pn Therefore, since p"(«i ) < p™ _1 (ui ), it follows by (72) that 



/r(«b) > , Lg, P"k) ^ ^ dcgl ( 73 ) 
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Since wo has been chosen arbitrarily, this completes the proof of the theorem. ■ 

The lower bound for the fixation probability in Theorem|4]is almost tight. Indeed, if a graph G = (V, E) 

with n vertices is regular, i.e. if degu = degu for every u, v € V, then f r (G) = ^_ I by Lemma JtJ (cf. also 

the Isothermal Theorem in [l6] ), and thus f r {G) = for large enough n. On the other hand, Theorem [4] 
implies for a regular graph G that f r (G) > £=p 



4.2 A class of selective suppressors 

In this section we present for every function 4>{n), where 4>(n) = w(l) and <fi(n) < ^/n, the class GM n ) — 
{G4,(n),n : n > 1} of , )-selective suppressors. We call these graphs 4>(n)-urchin graphs, since for 

4>{n) — 1 they coincide with the class of urchin graphs in Section 



3.2 



For every n, the graph G^ n \ n 
and V?, 



(y4>(n),m has n vertices. Its vertex set V^^ n y n can be partitioned into two sets V^ n ^ 

where \V}, 

[fif.it' (p{'H)-\-i 1 ip\iL),ib< (p^n)-\-L ' ipyiL),- 

has 4>{n) neighbors in Vh n ) n , and every 



(n),nl — <p{n) + l 

independent set in G0(„) jT1 . runuumuic, cvcij vhka u, c r M n ) t n llcia fVJ ^cignuuio m r^,^^ 
vertex v £ ^( n ) n has <fi 2 (n) neighbors in V^i n ) „■ Therefore degi> — n + 4> 2 {n) — 1 for every v € Vh n ) n ano - 



Furthermore, every vertex u G V? 



degu = 4>{n) for every w € K*(n) n- ^n example of a graph G0(„).„ is illustrated in Figure 10(a) 



0(n),n ■ 




• * " ' • V 2 , X 



.f^. -indep, set • 



7i ii 



^2 . 



^-2 * 



5*0 A <3l * S2 




(a) 



(b) 



Figure 10: (a) The graph G^i n \ n with n vertices and (b) the relaxed Markov chain. 



Lemma 8 For every v g „ sufficiently large n, 

fr(v) < 5, • ^ 
n 

Proof. Denote by the state, in which exactly > vertices of V^> n > n are infected and all vertices of 
^(n) n are n °t hifected. Note that Sq is the empty state. Furthermore senote by Ff. the state where exactly 
k > vertices of Vh n ) n ano - a ^ least one vertex of V^ n ) n are infected. In order to compute an upper bound 
for the fixation probability f r (Si) (i.e. of the fixation probability f r {v) where v £ Vh n ) «)' we can se ^ ^he 
value f r (S^) and the values / r (Ffc) for every k > 1 to their trivial upper bound 1. That is, we assume that 
the state S 1 ^, as well as all states F k , where k > 1, are absorbing. After performing these relaxations (and 
eliminating self loops) , we obtain a Markov chain, whose state graph is illustrated in Figure |10(b)| For any 
1 < k < 5 — 1 in this relaxed Markov chain, 

fr(S k ) = a k f r (S k+1 ) + Pkfr(Sk-l) + Ik (74) 
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where 



rk 



Oik 



Pk = 



Ik 



n—k 
n+0 2 (n) — 1 

<t> 2 ( n ) I n—k 

4>{n) ' n+0 2 (n)-l 



rk 



(n) — l 
J2k 



where T]j, = rk , ? 3 A , + 

■i— 'ft; n-fcWn) — 1 



2—k 



0(n) "r" n+0 2 (n) 



n:) + rt n+^w-r Note now b y < 75 ' that 



A = 1 ^ + 2 (») (n + ^ 2 (n)-l) \ > > : 

4>(n) (n — k) J r 



oik r 



Furthermore, since k < f , it follows by (75) that 



Tfc = cj) 2 (n) 2 4> 2 {n) 
oik n — k n 



Now, since + fik + 7fe = 1 and f r {Sk) > fr(Sk-i) f° r every fc, ([74]) implies by ( |76| and ([77]) that 

fr(Sk+l) - fr(Sk) 



— (fr(Sk) ~ /r(Sfc-l)) - "(I - /r(S*)) 
Oik Oik 



> 
> 

> 



— (/r(5fc) - /r(Sfc-l)) - 2^ 



(n) 



0(n) 



(/r(Sl)-/r(5o))"2^-X; 



/r(5l)-2 



0(") 



Thus, since f r (Sik) = 1 in the relaxed Markov chain, we have that 



Therefore 



and thus 



WrOSl) = X)(/r(5k+l)-/r(5*)) 



/r(Sl) E 



fc=l 
f-1 

> E 

^=1 



' At \\ k 

<P\n) 



<f>(n) 



k=0 



< 1 + 2- 



fr(Sl) - 2 



</>» 



</>» 



E 



- 1) £S 



0(n) 



- 1 



/r(Sl) < 2 



</> 2 (n) 



v n/2-l /0(n)y 



2r 



</> 2 (n) 



n(4>(ri)-r) n /2-i(^n)\ k 



2G 



Therefore, since 4>{n) = w(l) and r is constant by assumption, it follows that r < for sufficiently large n, 
and thus 

frfr) < 4r^M + I < &._*(») 



n n n 



Using Lemma [8] we can now prove the next theorem. 

Theorem 5 For every function <p(n) , where 4>{n) = w(l) and<p(n) < \fn, the class QM n ) — {G<t>(n),n : n > 1} 
of 4>(n) -urchin graphs is a class of ( » 4>{nj^~ se ^ ec ^ ve su PP ressors - 

Proof. It follows by Lemma 8 that, if u G VU •> n , then / r (u) < 5r^^ = A- , for any r > 1 and sufficiently 



n/4>(n) 

large n. Therefore, since |V£^ n | = 0(„" +1 for every graph G^( n ) jn , it follows by Definition that the class 
of graphs is a class of , ^y)-selective suppressors. ■ 
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